CMB Power Asymmetry from Non-Gaussian Modulation 
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Non-Gaussianity in the inflationary perturbations can couple observable scales to modes of much 
longer wavelength (even superhorizon), leaving as a signature a large-angle modulation of the ob- 
served cosmic microwave background (CMB) power spectrum. This provides an alternative origin 
for a power asymmetry which is otherwise often ascribed to a breaking of statistical isotropy. The 
non-Gaussian modulation effect can be significant even for typical ~ 1CP 5 perturbations, while 
respecting current constraints on non-Gaussianity, if the squeezed limit of the bispectrum is suffi- 
ciently infrared divergent. Just such a strongly infrared-divergent bispectrum has been claimed for 
inflation models with a non-Bunch-Davies initial state, for instance. Upper limits on the observed 
CMB power asymmetry place stringent constraints on the duration of inflation in such models. 



Large-scale features in the cosmic microwave back- 
ground (CMB) offer interesting avenues for testing phe- 
nomena that occurred at very early times in the Uni- 
verse's history. While most inflationary models pre- 
dict approximately scale- invariant, Gaussian fluctua- 
tions, some amount of non-Gaussianity is invariably gen- 
erated [J]. In this Letter, we show that even for an almost 
scale- independent power spectrum of curvature perturba- 
tions (i.e. <~ 1CU 5 in amplitude on all scales), primordial 
non-Gaussianity can lead to interesting, significant effects 
on the CMB, in particular a large angular scale modula- 
tion of the small-scale power spectrum. This is achieved 
without violating stringent observational bounds on non- 
Gaussianity in the subhorizon perturbations. 

There are some observational indications for a dipolar 
modulation of the CMB power spectrum [2-4]. Such an 
anisotropic CMB sky can be described by [5] 



6(n) = [I + /(n)]e(n), 



(1) 



where O(n) is the observed, anisotropic temperature fluc- 
tuation ST/T, while 0(n) is a statistically isotropic tem- 
perature field, and /(n) is the modulating function. Note 
that while 9(n) is statistically isotropic and is thus (sta- 
tistically) invariant under a rotation of the coordinate 
system, /(n) depends on fixed directions on the sky. 

The lowest order modulation is a dipole, as any 
monopole of /(n) is absorbed in the angle-averaged CMB 
power spectrum. The most recent analysis of [2] obtains 
a statistically significant dipolar asymmetry, while the 
WMAP team does not confirm this finding [6]. Hanson 
ct al. [7] find that beam asymmetries provide an expla- 
nation for the nonzero quadrupolar asymmetry. Several 
scenarios have been proposed in the literature to explain 
possible power asymmetries: [8, 9] considered remnants 
from the preinflationary phase, [10-12] proposed a single 
large-scale curvature perturbation, while [13] studied a 
spacelike vector field. These scenarios either involve a 
change in the inflation field Aip ~ A across the present 
horizon, many orders of magnitude larger than expected 
from the amplitude of fluctuations, or a breaking of the 



symmetries of the background. 

Alternatively, one can interpret a large-scale modu- 
lation of the CMB temperature fluctuations as due to 
a non-Gaussian coupling between long and short wave 
modes [14, 15]. The power spectrum of the Bardeen po- 
tential (j) on short scales is modulated by the presence of 
long modes if the fluctuations are non-Gaussian. We can 
Taylor expand the power spectrum of short modes (k) in 
the presence of long modes (ke <C k) [1] : 
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where P™ od (k) is the modulated power spectrum and 
Bff, is the bispectrum [33]. G(k, k^) can be understood 
as a scale- and orientation-dependent generalization of 
the dimensionless nonlinearity parameter /nl- 

The scenario we are considering is not statistically 
anisotropic in any fundamental sense; rather, the ob- 
served power spectrum P™ od (k) depends on the direction 
of k because the long modes in our particular realization 
of the Universe statistically pick out certain directions 
hi, and non-Gaussianity couples these long modes to the 
observable ones. Also, this effect doesn't require having 
a large amplitude long wave mode 0(k^); a large kernel 
G in the squeezed limit is sufficient. 

Observational bounds on primordial non-Gaussianity 
are rather tight [16], which might lead one to expect that 
the proposed effect must be small. The key point is that 
current observational constraints come from modes where 
both k^ and k are within our horizon. This is however 
not necessary for Eq. (2) to apply, allowing even super- 
horizon modes kg which we cannot directly measure to 
have an impact on observable modes k in the form of an 
anisotropic modulation. Two conditions should be met 
for this effect to be interesting: (1) the kernel G should 
be anisotropic, i.e. a nontrivial function of k • k £ ; and 
(2) G has to grow in the squeezed limit, i.e. scale like 
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k/kf to some positive power. Existing constraints effec- 
tively bound G only for moderate ratios of k/kg, while 
the superhorizon modulation effect is sensitive to larger 
ratios. We will interpret claims of power asymmetries in 
the literature as upper limits and use them to constrain 
models with such a strong coupling between short and 
long modes. 

Power asymmetry: The fluctuations of a statistically 
isotropic Gaussian field O(n) are specified through the 
spherical harmonic coefficients, (0z m 0;* m /) = 5u>6 mm >Ci. 
Adopting the notation of [17], the 0; m are related to the 
Bardeen potential perturbations 0(k) via 
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where A;(fc) is the photon temperature transfer function. 
The power spectrum of the Gaussian temperature fluc- 
tuations is then given by, 

C ( l ) =IJ k 2 dkP^k)\Mk)\ 2 . (4) 

On the other hand, the spherical harmonic coefficients of 
the modulated field [Eq. (1)] are 

0/m — 0(m = ^ 0;'m'/iM / ^^^taiilf^'ni'i 
LM,l',m' 

where we have expressed /(n) in terms of its multipole 
moments (with respect to a fixed coordinate system). 
The integral over three spherical harmonics can be writ- 
ten in terms of Wigner 3-j symbols, leading at linear 
order in f LM to 
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The 3-j symbols entail that l+l'+L is even, m' —m+M = 
0, and that \l — l'\ < L < I + I'. The latter condition is 
particularly relevant since we are interested in the case 
where L is much smaller than I, I'. Eq. (5) gives the 
covariance matrix of in multipole space in terms of the 
(fixed) multipole moments Jlm and the statistics of 0. 
As expected, the covariance is not diagonal, but it is very 
close to diagonal for 1,1' » L, i.e., it is nonzero only if 
< L. 

Non-Gaussianity: We assume that there is some gen- 
eral non-Gaussianity described to leading order by a 
bispectrum B^. We are interested in the limit kg -C 
Ho < k, where H is the Hubble scale today. Following 
Eq. (2), we expect that the presence of long- wavelength 
modes together with the mode coupling induced by non- 
Gaussianity lead to a breaking of statistical isotropy 



through the preferred direction k^. Consequently, we now 
calculate the covariance of the temperature field given 
Eq. (2). Multiplying Eq. (3), with 0;< m ', and integrating 
over one of the momenta leads to 
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where we set |k — k e \ ~ k in the squeezed- limit approxi- 
mation (corrections scale as kg/k and higher) . We obtain 

(0; m 0*/ TO ,) = Su'S mm iCt 
+ I ^T.^MCu^LMike), (7) 
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where we have defined 

C lv {ki) = \ J k 2 dk [A,(fc)A?,(fc) + A?(fc)A r (fc)] 
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using the fact that the kernel G only depends on the 
angle between k and k^. Comparing with Eq. (4), we 
see that Cu(ki) is equal to the temperature power spec- 
trum obtained when replacing P^(fc) — > Gl(/c, fc^P^fc), 
i.e. with a different initial power spectrum of curvature 
fluctuations. Thus, apart from the fact that the non- 
Gaussian covariance involves Cw, instead of Cu +Ci>i', 
it is identical in structure to the covariance obtained for 
the anisotropic field Eq. (5) [34]. The fractional differ- 
ence between Cw and Cu,Ci>i> is of order L/l <C 1. We 
will thus approximate Cw in Eq. (7) with (Cu +Ci'i')/2. 

We conclude that if Gl(A;, kg) is significant in the limit 
ki/k — > for some L > 0, the temperature fluctuations 
of the CMB appear as if they experience a (large-angle) 
modulation of multipole order L. In particular, this ne- 
cessitates an anisotropic coupling of long- and short- 
wavelength modes. We now calculate the amplitude of 
this modulation. For scale-free bispectrum shapes, the 
kernel moments in the squeezed limit (k t <C k) can be 
written as 



G L {k,k e ) = g L 



(10) 



where gL is a constant and oll gives the scaling in the 
squeezed limit. We also define the temperature power 
spectrum with a tilted spectral index n s — > n s + a, 
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where ko = 0.05 Mpc -1 is the pivot scale used to normal- 
ize Pj)(k). By comparing Eq. (7) with Eq. (5), we can 
then read off the relation between the long-wavelength 
perturbations and the anisotropy coefficients /lm, for a 
given / considered: 
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The multipole coefficients which give the amplitude and 
direction of the modulation are thus related to the given 
realization of the large- wavelength modes 4>Qce). The 
last factor in Eq. (12) encodes the fact that in general 
this modulation is Z-dependent; i.e. one is effectively 
adding a tilted CMB power spectrum C;(— at) with an- 
gular modulation to the angle-averaged CMB power spec- 
trum. While we cannot predict the direction of the power 
modulation, we can calculate the expectation value of the 
amplitude, defined as A = (^m=-l IIlm] 2 ) 1 ^ 2 - Since 
the Jlm are proportional to </>(1q), they are Gaussian- 
distributed complex numbers with mean zero. The am- 
plitude A thus follows a x distribution for 2L + 1 degrees 
of freedom, with an expectation value of 
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If at is sufficiently negative, (A) diverges as we let the 
lower integration bound go to zero. Such a prediction 
for A can be ruled out to high significance by the data 
if the observational limit An m <C (A). In general, if 
P<p(ki) oc fc" s ~ 4 , then at < (1 — n s )/2 for some L > in 
Eq. (10) is necessary for a significant large-scale asym- 
metry of the CMB. Fig. 1 shows quantitative results for 
the expected asymmetry (A) with L = 1, as a function 
of the CMB multipole I. We adopt the ansatz Eq. (10) 
with gi = 1, and integrate from fc^ min to fe«, max = l/rj\ ss , 
where i]i ss is the comoving distance to the last scattering 
surface (the latter choice is unimportant numerically). 
We choose three different sets of (a±, fc^min) and use 
CAMB [18] for the computation of C;(a). Clearly, a sig- 
nificant amplitude of power asymmetry can be achieved 
with a range of parameters. The closer a is to zero, the 
smaller fc^ m i n needs to be to generate a given amount 
of asymmetry (at fixed <?i). On the other hand, a more 
negative a.\ leads to a stronger scale dependence: the am- 
plitude of the modulation approximately scales as l~ ai . 

Inflationary bispectra which consist of symmetrized 
polynomials in the three momenta k\ , fo, &3 do not lead 
to a power asymmetry since the coupling of modes is 
isotropic (Gl = for L > 0). However, these simple 
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FIG. 1: Expected amplitude (A) of a dipole modulation (L = 
1) as a function of the CMB multipole I, for three sets of 
values for (ai, fc^min) as indicated in the figure. We have 
used Eq. (10) with gi = 1. Predictions for different L can be 
obtained by multiplying with g L (L\) 2 A L /2(2L)!. 



bispectra are often only obtained as separable approxi- 
mations to the more complicated exact bispectra, which 
may themselves in fact lead to Gt ^ 0. Hence, it is 
crucial to consider the full, exact bispectrum when de- 
termining whether a given inflationary model leads to a 
power asymmetry. It is clear however that such a power 
asymmetry requires a violation of the standard consis- 
tency relation [1], at least on the scales of interest, as it 
contains no anisotropic coupling between long and short 
modes. A recent example is solid inflation [19], which 
predicts a quadrupolar coupling between long and short 
modes. But since in this model G does not grow in the 
squeezed limit, the resulting quadrupolar modulation of 
the power spectrum is small. 

An example of a model that does produce a large-scale 
power modulation is the ekpyrotic ("case II") scenario 
of [20], which generates non-Gaussianities that in the 
squeezed limit lead to ot\ — — 1 — e and ai = — e, where 
e > is a red tilt. Thus, in this model one has diver- 
gent dipole and quadrupole modulations. Another case 
which has attracted recent interest is modifications to 
the initial state (non-Bunch-Davies) in single-field slow- 
roll inflation. These can lead to non-Gaussianity with 
a,\ = — 1 [21-27]. The squeezed bispectrum in the simple 
non-Bunch-Davies state considered in [26] reads 



£ (|k + k f /2|,|-k + k^/2|,|-kf|)=BP (^)P (fc) 



k( 



x Rc 



- 1 — p i(l+ju)fei/fe* _ 1 _ p i{l-n)k t /k„ 

/l^— + f2—, 

1 + /i 1 - jU 



(14) 



where Re(/i + /2)/2 ~ Nk, and N k is the occupation 
number of the momentum state k [35] , // = — k • k^, 
fc* ~ l/l^inl is related to the conformal time at which 
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the initial state is specified, and kg > fc* in order for 
this result to apply. B is a dimcnsionless constant equal 
to 4e in the case studied in [26, 27], although it could 
take larger values in more general models. The kernel 
Gl scales as k/kg in this model, with g^ oc Nk for even 
L and aj, = — 1. For odd L, the modulation scales as 
K~H which is suppressed by k e /k. We use observational 
upper limits on the (primordial) quadrupole modulation 
amplitude A < 0.1 [7] to place constraints on £^ im i n = 
Numerical evaluation of Eqs. (13)-(14) leads to a 95% 
C.L. lower limit of [36] 

fc # >2xl0" 5 /iMpc" 1 A r )c 6, (15) 

implying no more than <~ 3 e-folds of inflation beyond 
our current horizon for Nj-B ~ 1. This complements the 
bound on a non-Bunch-Davies initial state from backre- 
action arguments, which is sensitive to Nk but not B. 

Conclusions: Large-scale modulations of the CMB 
temperature fluctuations offer an interesting testing 
ground for the physics of the very early Universe. 
We have shown that certain types of primordial non- 
Gaussianity generically predict large power asymme- 
tries in the CMB. The requisite non-Gaussianity can be 
thought of as an anisotropic, scale-dependent /nl which 
grows in the squeezed limit. 

Upper limits on such a modulation can put strin- 
gent constraints on this class of models, which includes 
scenarios with a non-Bunch-Davies initial state. One 
can roughly estimate the modulation amplitude from 
the dimensionless bispectrum amplitude G(k, k^) for the 
longest observable mode kg <~ Hq 1 through 

(A) ~ G (k, h = H^ 1 ) 4 x 1(T 5 ( -^-) " , (16) 



where fc^ m i n refers to the longest superhorizon mode re- 
sponsible for the modulation, and controls how fast 
G grows in the squeezed limit (Eq. 10). Conversely, 
observational hints of a power asymmetry provide mo- 
tivation to further investigate such models. A convinc- 
ing detection of a CMB power asymmetry, if interpreted 
in terms of this scenario, would open an observational 
window to scales much larger than the present horizon 
(kg <C l/r7i ss ), which are otherwise completely inacces- 
sible to direct observation. This fact distinguishes this 
effect from a modulation of the temperature power spec- 
trum by a horizon-scale mode. 

We have shown that the power asymmetries are gener- 
ally scale dependent and increase toward smaller scales. 
Thus, unless one invokes a change in the shape or am- 
plitude of non-Gaussianity on smaller scales, a nonde- 
tection of a similar power asymmetry in the large-scale 
structure [28, 29] puts further stringent constraints on 
models that produce such asymmetries. Furthermore, 
models with bispectra that peak more strongly in the 
squeezed limit than the local model will in fact generate a 
scale-dependent bias in large-scale structure tracers [30- 
32] A6 oc k~ n with n > 2 [26, 27]. Observations of the 
large-scale structure will thus be of great importance in 
strengthening constraints on the possible non-Gaussian 
origins of a power asymmetry. 
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